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ABSTRACT 


This paper considers a four-bar linkage of which 
two adjacent hinges have skew axes and are fixed in the 
coordinate system. The remaining two joints are universal. 
Analytical relations are found between the rotations, 
angular velocities, angular accelerations, and torques 
about the fixed joints. 


I. INTRODUCTION 


In considering a spacecraft design having four solar panels arranged in two pairs, it was desired to 
investigate the possibility of utilizing a single, linear, spring actuator to unfold each pair from its launch 
position. Figure 1 shows a model of the spacecraft with one pair of panels folded against the spacecraft and 
one pair unfolded. The actuation mechanism consisted of a dual four-bar linkage, with skewed axes, attached 
to a single actuator. Figure 2 is a model of the linkage. The actuator is represented by the T-shaped handle 
on the right, and the solar panels by the section of beam and flat plate on the left. One panel has been dis- 
connected from the linkage for clarity. A linear motion of the actuator, then, transmits, through the linkage, 
a torque about each hinge axis H. 

Figure 3 diagrammatically depicts one-half of the linkage and defines the terms used in the analysis. 
Axes A and H are fixed to the spacecraft and are skewed relative to each other. Joints P and Q are' 
universal bearings. 
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II. DISCUSSION OF PROBLEM 


It was necessary to find the forces in the links and the torque relations between the arms. If the 
torques about A and H are both known, it is a simple matter to find the forces in the various links. Therefore, 
this paper concerns itself primarily with finding the relations between the torques about the two fixed axes. 

It is also intended that a dynamical analysis be facilitated, so the relations between the accelerations of the 
arms about the fixed axes are also determined. By-products of the investigation are the relations between the 
an gles and rates about the axes. 

Four-bar linkage solutions in the literature often consider the motion of points connected to link r. 
This paper does not, and the omission is part of the reason that the paper’s title refers to a “partial solution” 
to the problem. Since the locations of P and Q are known as functions of yS or S from the present results, 
the loci of such points may be found without unreasonable difficulty. 
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III. SOLUTION 

The torque relations will follow directly from the principles of Virtual Work and Conservation of 
Energy once the velocity relations between links c and d are known. 

The velocities and accelerations may be found by successive differentiations of the angles with 
respect to time. Consequently, it is first necessary to geometrically relate the angles /3 and 8. 

The location of point P can easily be found in terms of the Cartesian coordinates of Point A, the 
angle 8, and the length d. Similarly, the coordinates of point Q can be found from point H, the angle ft, and 
length c. The angles /3 and 8 may be combined in one equation by equating the distance between points P 
and Q to the length of link r (Eq. 1). 

After this equation is solved for /3 (Eq. 8) or 8 (Eq. 9), first and second time derivatives may be 
taken, giving velocity and acceleration relations. The torque calculation is inserted between the first and 
second differentiations, since it follows logically at that point. 

If the coordinates of A and H are given in some coordinate system such as in Fig. 3, say A % , A , 

A z and H x’ H z’ then 


b = H - A , 

X X x’ 


b = H - A , 

y y y’ 


b = H - A 

z z z 


The components of the location of point P, with a secondary coordinate system at A, are 


P % = 0, Py - -d cos 8, P z = -d sin 8 


The components of the location of point Q are 


Q x = b x - c sin a cos J3, + c cos a cos /3, Q z = b z + c sin /3 


The square of the distance between P and Q is 


<■ p , - V 2 * < p r - e r > 2 - - <?,) 2 ■ ' 


in 
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Expanding, Eq. (1) becomes 


(— b x + c sin a cos /3) 2 + (— d cos S ■ — — c cos a cos /3) 2 + (- d sin 8 — b z — c sin / 3 ) 2 = r 2 


(2) 


Further expansion gives 


6 2 — 26 c sin a cos B + c 2 sin 2 a cos 2 B + c / 2 cos'* 8 + 6 2 
xx ~ ~ y 

+ 2 by d cos S + c 2 cos 2 a cos 2 /3 + 2 cd cos a cos /3 cos S 

+ 2 be cos a cos B + d 2 sin 2 8 + 6 2 + 26, (/sin S + c 2 sin 2 B 
y ' z z ' 

+ 2 cd sin j8 sin 8 + 2b z c sin /3 = r 2 (3) 

Collecting terms in Eq. (3) gives 

a j cos /3 + ag si n /3 = 03 (4) 

or 


6 ^ cos S + 6 g s i Q S = 63 


where 

a, = 2c [-6 sin a + (6 4 - c? cos 5) cos a] 
lx y 

a 2 = 2c (6 z + (/ sin S) 

ag = r 2 - 6 2 - c 2 - (/ 2 - 2(/ (6 2 sin 8 + b^ cos 8) 
b j = 2,d{by + c cos a cos /3) 

&2 = 2c? (6 2 + c sin /3) 

6 g = r 2 - 6 2 - c 2 - d 2 - 2c [(6 y cos a - b x sin a) cos fi + b z sin /3] 
Equation (4) is considered. If 


A = 



+ a 


2 

2 


(5) 
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and 


then 


€ = tan 


a 

a 


1 

2 


A sin (/3 + e) = a 3 


(6) 



(7) 


a 3 o j 

J3 = sin -1 tan -1 (8) 

+ y/a\ + of a 2 

Equation (8) gives four, not necessarily distinct, possible solutions to Eq. (4). Physical considera- 
tions indicate that two, and only two, values will satisfy the original equation. Unless additional criteria 
are available by which one of the values may be eliminated, both must be carried throughout the remainder of 
the calculation. 

Since Eq. (5) is identical in form to Eq. (4), 




tan 


-1 


(9) 


2 ^ .2 

+ b 2 


Differentiation of Eqs. (8) and (9) gives the relations between the angular velocities. 


d/3 

dt 


2 2 

°1 + °2 


da„ 


dt 


da. 


dt 


A' 


da 0 


da. 


dt 


2 2 

°i + 4 


dt 


+ a. 


da n 


dt 


( 10 ) 


where 


a'= + 4 - 4 
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Similarly, 


dS 1 

L ib * 


± J_ 

db s h 

L ib > 

* 2 \ 
+ b 1 

dt b\ + b\ 

\ 1 dt 

°2 

dt y 

B' 

dt bj + &2 

V * 

+ °2 

dt ) 


( 11 ) 


whe 


yf l f 


2 + hi - bl 


Physical considerations indicate that A ' and B 'are real unless (1) the linkage is impossible, or (2) the 
oscillation range of the oscillating bar (if there is one) is exceeded. 

Since jd and 3 are the only variable quantities in the definitions of the quantities o and b, 

da \ d 8 

= —2 cd cos a sin 3 


dt 


dt 


(12) 


da 0 d 3 

_ = 2 cd cos 3 

dt dt 

da 3 d 3 

= 2d ( b sin 8 — b z cos 3) 

dt dt 


(13) 


(14) 


db. 


d ) 3 


= -2 cd cos a sin fd ■ 

dt dt 


(15) 


db 2 d/d 

= 2 cd cos fd 

dt dt 


(16) 


db i , dfd 

= 2c [(3 cos a — b sin a) sin fd — b cos /3J 

dt 7 dt 


(17) 


In order to find the torque relations at the hinges A and H, the fact that the work being put in at, say, 
hinge A equals the work taken out at hinge H (neglecting friction in the linkage) is used. If the torque at A 
is called T A , and the torque at H called T then, mathematically, this statement is 


6 
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T a d S = T h d/3 


(18) 


For generality, to avoid the problem of what happens if dS - d/3 = 0, dS and d/3 should be considered 
to be virtual displacements, and the works being equated are the virtual works. 

By virtue of Eqs. (12) -(14), it is seen that Eq. (10) may be written 



d S 

T 

“ 1 A/H 


dt 

dt 


(19) 


whe 


' A /H 


2 cd 


a l + a 2 


(< 2 j cos S + a 2 cos a sin S) 


2d 

A' 


o 

b sin S — b cos S — (a 0 cos S — a. cos a sin S) 

y 2 o o 1 t 

°1 + a 2 


(20) 


which depends entirely upon the geometry. Similarly, Eqs. (15) -(17) used in Eq. (11) imply 


d 8 
dt 


H / A 


d/3 

dt 


( 21 ) 


where 


2 cd 


1 H/A 


b\ + b\ 


(hj cos /3 + b 2 cos a sin j3). 


2c 

B‘ 


b 3 d 


( b cos a - 6 sin a) sin /3 - 6 cos (3 - (6, cos B - b, cos a sin B) 

J X Z _ r> O “ *■ 


b{ + b\ 


( 22 ) 


which depends entirely upon the geometry. Considering the differentials as virtual velocities, Eqs. (19) and 
(21) are statements that virtual power, in the absense of losses, is the same at both ends of the linkage. Then 


T A/H = torque ratio, torque at A divided by torque at H 
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and 


Tfl/A = torque ratio, torque at H divided by torque at A 


A second time differentiation of the expressions for /3 and S gives the relations between the angular 


accelerations of the arms. 



do 2 
dt 
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Regarding the ± sign, the same sign is used in Eqs. (23) or (24) as was used in Eqs. (10) or (11). 
Finally, the second derivatives of the a and b quantities and the derivatives of A 'and B 1 are as 

follows: 


dA ' 

1 

/ da j 


da 3 

— 

— 

(°1 

+ a 2 

a 3 

dt 

A' 

\ dt 

dt 

dt 


(25) 
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d 2 b l 
dt 2 


— 2 cd cos a- 


COS /3 



+ sin /3 



(30) 


dt 2 


- led 


sin /3 




d 2 b 3 
dt 2 


2c 



cos a — 


sin a) 


cos /3 





(31) 


(32) 
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NOMENCLATURE 


A 

a V a 2’ °3 

b v ^ 2 ’ b 3 
A, 6, A', B 

b, c, d, r 

H 


P, Q 

T „ 


l H 


a 


A § 


“driver” axis 

quantities defined for manipulative purposes 

lengths of links (Fig. 3) 

“follower” axis 

universal joints at ends of connecting rod 
torque about A 
torque about H 

angle between the axes A and H 

angles locating the positions of the “follower” and the “driver” (Fig. 3) 
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Fig. 1. Solar panel 
arrangement on 
spacecraft 
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note: no lengths or angles in the above views are 

TRUE UNLESS THEY ARE SO INDICATED. 


Fig. 3. Three-dimensional four-bar linkage 
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